Model Verification

Results: Mortality

ggplot(alltrtw1%>% filter(EP=="Mortality" & HT == "Fixed"),aes(x=Time,y=Mean))+annotate("rect", xmin = 5*30, xmax = 6*30, ymin = -Inf, ymax = Inf,alpha = .1,fill = "blue")+geom_line(aes(x=Time,y=Mean,col=factor(HCx)))+geom_point(aes(x=Time,y=Mean,col=factor(HCx)))+geom_ribbon(aes(ymin=`5th Per.`,ymax=`95th Per.`,fill=factor(HCx)),alpha=0.1)+facet_grid(HT~.)+ylab("Total Mass")+xlim(c(120,200))+ggtitle("Looking at Mortality Alone")+labs(color="Sensitivity Level",fill="Sensitivity Level")

Figure 1: Total Mass Response with Mortality Effect

Results: Effect parameters

ggplot(alltrtw1%>% filter(HT == "Fixed"),aes(x=Time,y=Mean))+annotate("rect", xmin = 5*30, xmax = 6*30, ymin = -Inf, ymax = Inf,alpha = .1,fill = "blue")+geom_line(aes(x=Time,y=Mean,col=factor(EP)))+geom_point(aes(x=Time,y=Mean,col=factor(EP)))+geom_ribbon(aes(ymin=`5th Per.`,ymax=`95th Per.`,fill=factor(EP)),alpha=0.1)+facet_grid(HCx~.)+ylab("Total Mass")+xlim(c(120,200))+ggtitle("Sensitivity: 0.05 vs. 0.95")

Figure 2: Total Mass Response with different Effects

Looking at each EP

ggplot(alltrtw1%>% filter(HT == "Fixed"),aes(x=Time,y=Mean))+annotate("rect", xmin = 5*30, xmax = 6*30, ymin = -Inf, ymax = Inf,alpha = .1,fill = "blue")+geom_line(aes(x=Time,y=Mean,col=factor(HCx)))+geom_point(aes(x=Time,y=Mean,col=factor(HCx)))+geom_ribbon(aes(ymin=`5th Per.`,ymax=`95th Per.`,fill=factor(HCx)),alpha=0.1)+facet_wrap(EP~.)+ylab("Total Mass")+xlim(c(120,200))+ggtitle("Sensitivity==0.05 vs. 0.95")+labs(color="Sensitivity Level",fill="Sensitivity Level")

Figure 3: Total Mass Response with different Effects

Above Mass

ggplot(alltrtw1%>% filter(HT == "Fixed"),aes(x=Time,y=Mean))+annotate("rect", xmin = 5*30, xmax = 6*30, ymin = -Inf, ymax = Inf,alpha = .1,fill = "blue")+geom_line(aes(x=Time,y=Mean,col=factor(HCx)))+geom_point(aes(x=Time,y=Mean,col=factor(HCx)))+geom_ribbon(aes(ymin=`5th Per.`,ymax=`95th Per.`,fill=factor(HCx)),alpha=0.1)+facet_wrap(EP~.)+ylab("Above Mass")+xlim(c(120,200))+ggtitle("Sensitivity==0.05 vs. 0.95")+labs(color="Sensitivity Level",fill="Sensitivity Level")

Figure 4: Total Mass Response with different Effects

Below Mass

ggplot(alltrtw1%>% filter(HT == "Fixed"),aes(x=Time,y=Mean))+annotate("rect", xmin = 5*30, xmax = 6*30, ymin = -Inf, ymax = Inf,alpha = .1,fill = "blue")+geom_line(aes(x=Time,y=Mean,col=factor(HCx)))+geom_point(aes(x=Time,y=Mean,col=factor(HCx)))+geom_ribbon(aes(ymin=`5th Per.`,ymax=`95th Per.`,fill=factor(HCx)),alpha=0.1)+facet_wrap(EP~.)+ylab("Below Mass")+xlim(c(120,200))+ggtitle("Sensitivity: 0.05 vs. 0.95")+labs(color="Sensitivity Level",fill="Sensitivity Level")

Figure 5: Above Mass Response with different Effects

Shannon

ggplot(alltrtw1%>% filter(HT == "Fixed"),aes(x=Time,y=Mean))+annotate("rect", xmin = 5*30, xmax = 6*30, ymin = -Inf, ymax = Inf,alpha = .1,fill = "blue")+geom_line(aes(x=Time,y=Mean,col=factor(HCx)))+geom_point(aes(x=Time,y=Mean,col=factor(HCx)))+geom_ribbon(aes(ymin=`5th Per.`,ymax=`95th Per.`,fill=factor(HCx)),alpha=0.1)+facet_wrap(EP~.)+ylab("Shannon")+xlim(c(120,200))+ggtitle("Sensitivity==0.05 vs. 0.95")+labs(color="Sensitivity Level",fill="Sensitivity Level")

Figure 6: Shannon Response with different Effects

NPFT

ggplot(alltrtw1%>% filter(HT == "Fixed"),aes(x=Time,y=Mean))+annotate("rect", xmin = 5*30, xmax = 6*30, ymin = -Inf, ymax = Inf,alpha = .1,fill = "blue")+geom_line(aes(x=Time,y=Mean,col=factor(HCx)))+geom_point(aes(x=Time,y=Mean,col=factor(HCx)))+geom_ribbon(aes(ymin=`5th Per.`,ymax=`95th Per.`,fill=factor(HCx)),alpha=0.1)+facet_wrap(EP~.)+ylab("NPFT")+xlim(c(120,200))+ggtitle("Sensitivity: 0.05 vs. 0.95")+labs(color="Sensitivity Level",fill="Sensitivity Level")

Figure 7: NPFT Response with different Effects